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Abstract 

In this paper, we consider general two-way relay networks (TWRNs) with two source and N relay 
nodes. A distributed differential space time coding with analog network coding (DDSTC-ANC) scheme 
is proposed. A simple blind estimation and a differential signal detector are developed to recover the 
desired signal at each source. The pairwise error probability (PEP) and block error rate (BLER) of 
the DDSTC-ANC scheme are analyzed. Exact and simplified PEP expressions are derived. To improve 
the system performance, the optimum power allocation (OPA) between the source and relay nodes 
is determined based on the simplified PEP expression. The analytical results are verified through 
simulations. 
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I. Introduction 

It is well known that cooperative communication improves system robustness and capacity by 
allowing nodes to cooperate in their transmission to form a virtual antenna array (T). Compared 
to one-way relay networks (OWRN), two-way communication is an effective scheme to improve 
the spectral efficiency by allowing the simultaneous exchange of two-way information flows. 
In 0, the authors first studied the two-way relay networks (TWRN) and derived its achievable 
bidirectional rate. The TWRNs have attracted increased interest due to its high spectral efficiency. 
Various protocols for the TWRNs have been proposed recently 01, [|4]|. 

In BH, the conventional network coding scheme was applied to the TWRNs. In this scheme, 
two source nodes transmit signals to the relay, separately. The relay decodes the received signals, 
performs binary network coding, and broadcasts network coded symbols back to both source 
nodes. However, this scheme may cause irreducible error floor due to the detection errors which 
occur at the relay node. In [3 J, an amplify and forward based network coding scheme was 
proposed. In this scheme, both source nodes transmit at the same time so that the relay receives 
a superimposed signal. The relay amplifies the received signal, and broadcasts it to both source 
nodes. Each source node subtracts its own contribution and estimates the signal transmitted from 
the other source node. Analog network coding is particularly useful in wireless networks as the 
wireless channel acts as a natural implementation of network coding by summing the wireless 
signals over the air. 

Recently, distributed space-time coding for OWRNs was proposed in to achieve spatial 
diversity. Since OWRNs take place only in a single-direction, to further improve the spectral 
efficiency of the relay networks, the distributed space-time coding was proposed for TWRNs 
in (H, 0. However, most of the existing works on distributed space-time coding in TWRNs 
consider coherent detection at each receiver with the assumption of available channel-state 
information (CSI). In some situations, e.g., the fast fading environment, the acquisition of 
accurate CSI presents great challenge, and training becomes expensive and inefficient while there 
are a large number of relays in the wireless networks [8J. hi this case, differential modulation 
would be a practical solution because it requires no knowledge of the CSI. 

The distributed differential space-time coding was first proposed for OWRNs in BH. In TWRNs, 
the signal received at the relay node is a superposition of two symbols sent from two source 
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nodes. Thus if there is no CSI available at source and relay nodes, it will be very difficult to 
design distributed differential modulation schemes in TWRNs. The challenge is due to the blind 
channel estimation from the superimposed signals at the relay and unknown self-interference at 
each destination. In [10], the authors first extended the distributed differential space-time coding 
to TWRNs. In order to enable differential encoding and decoding, this scheme starts with a 
four-stage initialization phase, which is similar to traditional one-way relaying, to transmit the 
bi-directional reference signals respectively. After initialization, each user then proceeds to the 
data transmission. Information exchange between two users is done in two time slots. However, 
the decoding algorithm in ifTOl is a non-coherent detection scheme where the decoding of current 
symbol is based on the estimation of the previous symbol. Consequently, when one symbol 
was decoded incorrectly, it will affect the decoding of consecutive symbols thus leading to 
serious error propagation. To solve this problem, periodical initialization of the protocol has 
to be performed to transmit new reference signals for decoding, making the proposed scheme 
inefficient. Furthermore, no pairwise error probability (PEP) analysis was performed in IfTOl 
due to the complexity of the protocol. [j8] presented an analog network coding scheme with 
differential modulation using the amplify-and-forward protocol for bidirectional relay networks. 
However, this scheme is limited to single relay node, thus cannot be extended to the distributed 
space-time codes. 

Unlike lf8l- lfT0"l . in this paper, we propose a distributed differential space time coding with 
analog network coding (DDSTC-ANC) scheme for the TWRNs with multiple relays. In this 
scheme, two source nodes perform differential modulation, and transmit the differential modu- 
lated symbols to all the relay nodes in the first time slot. The signal received at the relay node is a 
superposition of two transmitted symbols. In the second time slot, the N relay nodes broadcast the 
processed signals to both source nodes simultaneously. We propose a blind estimation technique 
which can be used to subtract the self-interference without knowledge of CSI at both relay 
nodes and two source nodes. A simple differential signal detector is then developed to recover the 
desired signal at each source. The performance of the proposed differential DDSTC-ANC scheme 
is analyzed and the PEP and block error rate (BLER) expressions are derived. They show that 
the proposed differential scheme can achieve the same diversity order as the coherent detection 
scheme but is about 3dB away compared to the coherent detection scheme due to the differential 
transmission. To further improve the system performance, the optimum power allocation (OPA) 
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between the source nodes and the relay nodes is determined based on the provided simplified 
PEP expression. The analytical results are verified through simulations. Simulation results also 
show that the proposed differential scheme with OPA yields superior performance improvement 
over an equal power allocation (EPA) scheme. 

The rest of this paper is organized as follows: In Section HH the system model is introduced. 
Section [Till presents the proposed DDSTC-ANC scheme. The performance and diversity order 
of DDSTC-ANC are analyzed in Section El In Section |Vl the OPA for the DDSTC-ANC is 
presented. Simulation results are provided in Section |VTJ In Section IVIIL we draw the main 
conclusions. 

Notation: Matrices and vectors are denoted using capital letters and boldface lowercase 
letters, respectively. (■)*, (-) T and {-) H represent conjugate, transpose and conjugate transpose, 
respectively, for both matrix and vector. For a complex matrix A, det A and ixA denote the 
determinant and trace of A, respectively. I m is the m x m identity matrix, diagjai, • ■ ■ ,a n } 
stands for an n x n diagonal matrix whose 2th diagonal entry is <Zj. In represents the natural 
logarithm, and || • || is the Frobenius norm. E and P(-) denote the expectation and probability, 
respectively. 

II. System Model 

In this paper, we consider a general TWRN with N + 2 nodes, as shown in Fig. [Q where two 
source nodes, T\ and T 2 , want to exchange information with each other through N relay nodes. 
It is assumed that each node in the network is equipped with one single antenna working in the 
half-duplex mode. We consider a quasi-static fading channel, where the channel remains constant 
for the duration of a frame and varies independently from one frame to another. Let fa and g, t 
denote the complex fading channel coefficients of 7\ — Ri and T 2 — Ri, respectively. Furthermore, 
we assume Rayleigh flat fading channels, i.e., /, ~ £/V(0, aj.) and g; L ~ £A/"(0, <t^), respectively. 
For analysis tractability, symmetry of the relay nodes is assumed in this paper, i.e., a - /. = (Tf, Vz 
and a 9i = cr g ,\fi. 

A general two-time slot TWRN protocol is used, as shown in Fig. \T\ In the first time slot, both 
T\ and T 2 transmit their messages and the relays {Ri, • • • , Rn} receive a superposition of the 
signals transmitted from T\ and T 2 . Let s(t) = [si(t), ■ ■ ■ , s T (t)] T and d(t) = [di(t), ■ ■ ■ , d T (t)] T 
denote the transmitted symbol vectors of Ti and T 2 at time t, respectively. They are normalized 
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as E{s(£)s(£) } = E{d(i)d(i) } = I T . The received signal vector at Ri can be written as 



n(t) = VPJi(t) s (t) + VP 2 9i(t) d(t) + Vi(t) 



(i) 



where Pi and P 2 denote the transmit power of T\ and T 2 , respectively, and Vj(t) represents the 
noise vector at Ri and each noise term follows a zero-mean complex additive white Gaussian 
distribution, i.e., v, ; (t) ~ CAf(Q,N I T ). 

During the second time slot, Ri processes Ti(t) to generate a space time coded symbol vector 
Xj(£). In this paper, we consider the amplify-and-forward protocol in the relay nodes. The transmit 
signal at the i-th relay is designed to be a linear function of its received signal and its conjugate 
CD: 



Xi(t) = Pi{t) (AiTi(t) + BiTi(ty 



(2) 



where A4 and Bi are two T x T complex matrices specifically designed for the construction of 
distributed space-time codings, and 0i(t) is the scaling factor at R4. 

In this work, the scaling factor /3i(t) in Eq. (0 can be obtained based on the available statistical 
CSI, which is specifically given by lfT2l 



Ri 



(3) 



where P R . is the transmitted power of R { . Since we assume P R . = P R , we have j3i(t) = /3,Vz. 
For simplicity, in this paper, we only design the system that either A4 is unitary, Bi = Ot (case 
I) or Bi is unitary, Ai = T (case II). Thus, case I means that the i-th column of the code matrix 
(S(t) and D{t) in Eq. ©) contains only the transmitted symbols, and case II means that the i-th 
column of the code matrix contains the linear combinations of the conjugate of the transmitted 
symbols only. Further more, we assume that T = N, i.e., the number of symbols in a space-time 
block code is equal to the number of relay nodes. We further define 



Oi = A i} fi = fi, g { = g h Vi(t) = V;(t), 
Si{t)±s(t), di(t)±d(t), ifBi = T , 

O^B h f { ±f*, g^gl Vj(t) = Vj(t)*, 
Si(t)^s(t)*, di(t) = d(t)\ if Ai = T . 



(4) 
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Then the relay node R { broadcasts the coded symbol vector Xj(£) back to both source nodes. 
Since Ti and T 2 are mathematically symmetrical, for simplicity, in the following, we only discuss 
the decoding and the analysis for the signals received by T 2 [fT3l . The received signal vectors 
at T 2 is given by 

N 

y 2 (t) = ^2gi(t)xi(t) + w 2 (t), (5) 

where w 2 (i) denotes the independent and identically distributed (i.i.d) additive white Gaussian 
noise (AWGN) vectors at T 2 and we have w 2 (t) ~ CJ\f(0,N Q l T ). 
The received signal at T 2 can then be rewritten as: 

y 2 (t) = V^5(t)h 12 (t) + y/F 2 D(t)h 22 (t) + n 2 (t), (6) 

where 

S(t) = [Oisi (*),••■ ,OjvSiv(t)], 
D(t) = [0 1 d 1 (t),---,0 Jv d JV (t)], 
< M*) = ■ ■ ■ J N (t)f N (t)g N (t)} T , (7) 

h 22 (t) = \Pi{t) gi (t) gi {t),--- ,(3 N {t)g N (t)g N {t)} T , 
= Eii ft(*)OiVi(t) + w 2 (t). 
It is easy to prove that E{n 2 (t)n 2 (t) H } = a 2 n2 (t)I N , and a 2 n2 (t) = (J2? =1 |A(t)| 2 |^(t)| 2 + l)iVo. 

III. Distributed Differential Space-Time Coding For TWRNs 

In this section, we propose a distributed differential scheme. First, we blindly estimate channel 
h 22 (t) defined in Eq. ©, which can be used to subtract the self-interference. Then, a simple 
differential signal detector is developed to recover the desired signal at source T 2 . 

In the proposed DDSTC-ANC, Ti encodes a message at time t into aniVxA^ unitary matrix 
U(t), which is then differentially encoded as s(£) = U(t) ■ s(t — 1), where s(t — 1) is the signal 
transmitted by Ti at time t — 1. Similarly, T 2 differentially encodes a message at time t into an 
N x N unitary matrix V(t), which is then differentially encoded as d(t) = V(t) ■ d(t — 1). 

For the first block, we can transmit a known vector to both source nodes that satisfies 
E{s(t) H s(i)} = E{d(t) H d(t)} = N, for example, [11 ■■• if or [VN ••■ 0] T . Similar 
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(8) 



to the differential space-time coding for multiple- antenna systems, having U(t) and V(t) unitary 
preserves the transmit power. 

For simplicity, we define Ui{t) = U(t) if B { = T , and U { (t) = U(t)* if A { = T . In the 
distributed differential scheme, the codes U(t) and V(t) should commute with the relay matrices 
[U, i.e.,Q 

oA(t) = u(t)Oi, 

or equivalently, 

AiU(t) = U(t)Ai, if B t = T , 
B i U*(t) = U(t)B i , if Ai = T . 

Hence, S(t) can be rewritten as 

S(t) = [0i#i(t)si(t - 1), ■ ■ • , N U N (t)s N (t - 1 

= U(t) ■ (OMt - 1), • • • , N s N (t - 1)) (9) 
= U(t) ■ S{t-1). 
Similarly, we have D(t) = V(t) • D(t - 1). 

The distributed differential space-time codes (STC) for TWRNs should be designed to satisfy 
Eq. ®. The design and choice of appropriate codes is beyond the scope of this work, here, 
we only briefly introduce some existing STCs that can be used in TWRNs. For the TWRNs 
with two relays, we can use Alamouti code lfT8l . which has full diversity and linear decoding 
complexity. Square real orthogonal codes (SORC), which have full diversity and linear decoding 
complexity, were proposed in [|9]| for two, four and eight antennas systems. 

Theorem 1: If the relay matrices have the property: tr{ OjOf} = N for i = j, tr{OjOf} = 
for i j, we have 

E{D(t) H y 2 (t)} = v^iVh 22 (t), (10) 
and h 2 2(t) can be approximated as 

1 1 L 

h 22 (t) « £ D( f - /)"y 2 (t - /), (ID 

^ 2 i=i 



More properties about the differential space-time coding can be found in |14|-|17|. 



s 



where L denotes the number of STC symbols in a frame. 

Proof: It can be proved by direct matrix multiplication and expectation. Due to the limited 
space, we omit the details. ■ 

We note that since receiver T 2 knows the symbols d(t) sent by itself, using the blindly 
estimated channel h 22 (t), we can subtract the self-interference at T 2 without using pilot symbols 
at the beginning. Although we can blindly estimate channel h 22 (t), T 2 doesn't have any CSI of 
hi 2 (£). Then based on the above theorem, a simple differential signal detector is developed to 
recover the desired signal s(t) at source T 2 . In the later performance analysis section, we assume 
that h 22 (t) is perfectly cancelled. Most of papers on distributed STCs for TWRNs also assume 
perfect self-interference cancellation, such as 0, lfl9l for coherent systems and (HI, lfl"3ll for 
differential systems. However, in practice, the estimation error will introduce some performance 
degradation which depends on estimation accuracy of hi 2 (t). The estimated h 22 (t) is used in 
simulations in this paper. In the simulation section, we have simulated the proposed scheme 
using the estimated h 22 (t) and the results show that the performance loss due to the h 22 (t) 
estimation error is negligible. 

By using Eq. ©, Eq. (fTTI) and the assumption of h 12 (t) = h 12 (t — 1), we have 

Y2(t)=y 2 (t)-^F 2 D(t)h 22 (t) 

= V^>(f)Mt) + n a (f) (12) 

= [/(t)y 2 (t-l)+n 2 (i), 
where n 2 (t) = n 2 (t) - U(t)n 2 (t - 1). Note that E{U(t)U(t) H } = I N , and n 2 (t) and n 2 (t - 1) 
are independent complex Gaussian random vectors with zero mean and covariance cr^ (t). We 
have E{h 2 (t)n 2 (t) H } = *i(*)Ir, where a\ 2 {t) = 2(£Zi Mt)\ 2 \9i(t)\ 2 + 1)N . Thus, n 2 (i) is 
a Gaussian random vector with zero mean and covariance <j\ 2 {t). 

Hence, the least square (LS) decoder can be performed to recover the transmitted signal 

arg min ||y 2 (t) - U k {t)f 2 (t - 1) || . (13) 

IV. Pairwise Error Probability and Block Error Rate Analysis 

In this section, we derive the PEP and the BLER of the proposed DDSTC-ANC scheme. 
Asymptotic diversity order is also analyzed in this section. 
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A. Pairwise Error Probability 

For simplicity, we define U^,kj{t) — U k (i) — Uj(t) and S^j(t) = S k (t) — Sj(t). The PEP 
of mistaking the fc-th STC block by the j-th STC block can be evaluated by averaging the 
conditional PEP over the channel statistics, i.e., fi,gi, and we havej [20J 



^(7)=E /4l * 




(14) 



2<W 

where 7 = P/N is signal-to-noise ratio (SNR), P is the total power in the TWRN and Q(x) = 
l/y/2n exp(— t 2 /2)dt is the Gaussian Q-function. Since it is very difficult to analyse y2(t — 1) 
directly, we approximate it using Eq. (fl~2l) as ^(t) ~ VP~l S(t)hi 2 (t). This approximation is 
particularly accurate at high SNR. Then based on Eq. ©, we have SA,ij(t) = UA,ij(t)S(t — 1) 
. We further assume hi 2 (t — 1) m hi 2 (t) . Then Eq. (fT4l) can be further simplified as 



Similarly, the PEP for the coherent scheme can be derived as 



Since er~ 2 (t) = 2cr^ 2 (t), the distributed differential scheme in TWRN is supposed to have 3dB 
loss in coding gain compared to distributed coherent scheme. 

Before deriving the PEP, we first define h 12 (t) = pG(t)f(t) , where f(t) = [fi{t), ■■■ , f N (t)] T 
and G(t) = diag{gi(t), ■ ■ ■ ,g N (t)}. Then we have the following lemmas. 

Lemma 2: The probability density function (PDF) of f(i) can be derived as 

p H = ^f^{- A ^ n )- (17) 

Proof: Since ~ CN(0,a}), we can prove that /?(*) ~ CAf(0,a 2 f ). Hence, ~ 
CjV(0,(T^). Note that /i(t), ■ ■ ■ , /jv(t) are independent, we can easily derive Eq. (fTTT) . ■ 
Lemma 3: B represents annxn Hermitian matrix (i.e., £? H = B), and x is an n x 1 complex 
vector. We have 



/ exp (-x H Bx) dx = Tr^der^S). 
jo 



(18) 



2 The superscript "d" denotes differential scheme and "c" represents coherent scheme. 
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Proof: Please see (21). ■ 
Note that the canonical representation of Gaussian Q-function is in the form of a semi-infinite 
integral, which makes analysis very difficult. Here, we use an alternative representation of the 



Gaussian Q-function from ll22l Eq. (4.2)] as Q(x) = ^ JJ^exp ^- 
some manipulations, we have 



2sin 2 e 



d9. Then, by doing 



Jtf,-(7)=E 



fi,g 



o 



' 71 
tt/2 

tt/2 



tt/2 



cxp 



2 sin 2 9 
E 9i [det(I + K'(9,t)]' 1 d9 



dtf 



(19) 



E, 



N 



n(i+^,*)AiiMor) 



i=l 



d0, 



where 1T(£) 



ft |/3| 2 G(t) H 5 A , fcj (f ) g 5A, fcj (t)G(t) 



,K'(9,t) 



a 2 f K(t) 



,l(0,t) 



Pi\ 



4(Ef=i |j8| a l»(*)l a +l)iVo 2 fl m"e' •V'"-' 8(Ef=i |^| 2 |fli(t)| 2 +l)iVo S in 2 e 

and Aj, z e {1, • • ■ , TV}, denotes the singular value of SA,kj(t) H S^kj(t)- The second step of the 
equation is based on the Lemma [2] and Lemma [3j 

Note that the mean of \gi(t)\ 2 is cr 2 . It is reasonable to approximate the term J2iLi Ifi'iWI 2 m 
1(9, t), by \9i(t)\ 2 ~ MtJ, especially for large JV (by the law of large numbers) 0, |fT2l 

Hence, 



l(6,t) w Z'(0) 



8(iV|/5| 2 (j2 + l)Ar sm 2 ^ 



(20) 



Let | ft (*)p 



^exp 



7i(*) 



7i(t). Since ^ ~ £A/"(0, cr 2 ), the PDF of ^(t) can be obtained as p(7»(t)) = 
. Hence, after doing some manipulations, the MGF-based PEP expression is 



derived as 



1 /"r/2 



7T 



E„ 



tt/2 AT p 

n 

i=l 



' N 

H{i + ne)x l \g l (t)\ 2 ) 

.1=1 
,sin 2 0, 



d0 



(21) 



where M 



Pl|/3| 2 ^a 2 
8(N\/3\ 2 aj+l)N 



i=i *- 

— — A,- and Ei(x) 



exp 

= J! 



sin 2 



Ei 



sin 2 

M 



d0. 



oo t 



dt, for a; < 0, is the exponential integral 



function ^ 8.211.1]. 

Next let us derive the simplified PEP expression at high SNR. Note that J23l 8.214.1] Ei (x) = 
C + ln(-ar) + J2T=i m> i x < °L where C is Eulers constant and C « 0.577 L23 , 9.73]. When 
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x tends to 0, the exponential integral function can be approximated as Ei (x) ~ ln(— x), for 



x < 0. At high SNR, we have exp ^ 
integral function, we have 



sin 2 6> 



1 and using the approximation for the exponential 



7T 



tt/2 N 

n 



,sin 2 ^ 



x 2 In (sin 9) + In 



(22) 



d6. 



" Mi ' \ \M„ 

i=i <- 

Note that £ /2 In sinxdx = — | In 2 ll23l 4.224.3]. Hence the ln(sin#) can be ignored, especially 
at high SNR. Using J23l 3.621.3], we have £ /2 sin 2 " 1 xdx 
simplified as 



(2 (2m)M !! f- The PEP can be further 



d l(2iV-l)H TT 



N 



(2N)l\ 



i=l 



' — ) In (Mi' 



(23) 



Finally, we derive the well-known Chernoff-bound-based PEP expression. From Eq. (fT9l) , 
setting 9 = tt/2, and doing some manipulations, the Chernoff-bound-based PEP expression is 
given as 



f(t) B K(t)f(t) 



1 N 

III 



(24) 



1=1 



' — )ln(Mi' 



The average BLER can be obtained based on the well-known union bound as 



(25) 



B. Diversity Order 

In this subsection, we analyze the asymptotic diversity order of the proposed DDSTC-ANC 
scheme. Firstly, we define the total transmission power is N ■ P. Note that N ■ P = N ■ P 1 + N ■ 
P 2 + N 2 ■ P Ri , P 1 = a x P, and P 2 = a 2 P. Denote the SNR 7 = P/N . Then, we rewrite M t at 
high SNR as M { = CXi 7, where C 



U-°l-°2) 2 2 
Qi N a S a 3 



(l-ai-a 2 ) 2 2 



8((l-oi-02)cr|+oio- 2 +a20"2+l/7) 8((l-«i-a2)cr|+«io-J+a2°J) ' 



Thus, the simplified PEP at high SNR can be rewritten as 



1 (2N 



1 



2 (2iV)!! Uf =1 CX 



n(ln(CA,)+ln( 7 )) 



i=i 



(26) 



(C 7)^^(7)] 
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i 

where C = \^ ^^r -^jv 1 cx j N ■ When S^kj(t)SA,kj(t) H is full rank, the diversity can be 
obtained as ll24l 

^h-^y^fi-^). (27) 

7^00 log(7) V log(7) J 

Thus, the diversity of the proposed DDSTC-ANC scheme for TWRNs is N(l— log log(7) / log(7)). 

V. Optimum Power Allocation 

In this section, we derive the OPA between the source nodes and the relay nodes that minimizes 
the total PEP in the TWRNs. Because the MGF-based PEP expression is very hard to analyze 
and gives little insight, we use the simplified PEP expression to derive the OPA. Here, we 
consider the total PEP in the TWRNs, and denote the PEP in T\ and T 2 as Pj 1 ^) and Pf- 2 {^), 
respectively. C in Subsection llV-Bl is rewritten as Gt x and Ct 2 for T\ and T 2 , respectively. Hence, 
we have 



13 V " " -« v "-2 (2N)\\ YH lXi (2g) 
x(C7^ + C7^)7- JV [ln(7)] JV , 

a2 (l-"l-"2) cr 2 (T 2 ai (1-^1-^2) ^2^ 

where CV ~ -^m rV; — V? and CV, ~ rV; — — ar- It is obvious that to 

minimize the PEP at high SNR, we should minimize the C£ + C^ 2 N in Eq. (|28T) .i.e., 

min^ + C^}, s.t.<; ^ (29) 

ai,«2 ' 

As a special case, when a 2 = a 2 = a 2 , we have oil 

Cr =C T = 2a{l - < — (30) 

Tl Ta 16JV " 647V' 

with equality when a = ~, or equivalently, P x = P 2 = j and Pr, = Thus, the OPA is 

such that the source nodes use half the total power and the relay nodes share the other half. We 

should emphasize that this power allocation only works for the TWRNs, in which all channels 

are assumed to be i.i.d. Rayleigh and no path-loss is considered. It is obvious that it may not 

be optimal when the path-loss effect is considered in the TWRNs. 

As the expression in Eq. (1291 is complicated, it is difficult to derive the closed-form solution 

for OPA when a 2 7^ a 2 . Here, we use numerical method, such as the nonlinear optimization 
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method, to obtain the optimal solution. In Section [V]] it is interesting to find that when a'j ^ cr|, 
cti + a 2 = 0.5 still holds for the simulated scenarios, which means the source nodes still share 
half the total power. 

VI. Simulations 

In this section, we provide simulation results for the proposed DDSTC-ANC scheme. Simu- 
lations are performed with PSK modulation and a frame size of 100 symbols over a quasi- static 
Rayleigh fading channels without specific mention. The estimated hn(i) and h 22 (t) are used 
in simulations. For comparison, we also present simulations over a GSM channel model with 
a symbol sampling period of T s = 3.693/is and a maximum Doppler shift of 75 Hz ifTOl . This 
ensures a slowly changing channel and allows the assumption of a constant channel over two 
consecutive time blocks. Without specific mention, we assume that cr? = a 1 = 1 and the source 
nodes uses half the total power and the relay nodes share the other half, i.e., Pi = \P, P 2 = \P 
and P Ri = ±P. 

From Fig. |2l we present the simulated BLER performance for the proposed DDSTC-ANC 
schemes using Alamouti for TWRNs. The performance of the corresponding coherent detection 
is plotted as well for better comparison. It shows that the differential scheme suffers about 3-dB 
performance loss compared with the corresponding coherent scheme, which has been validated 
in Subsection HV-Al Fig. |2| also compares the simulated BLER performance for our proposed 
DDSTC-ANC and the differential scheme in iflOll . It can be observed that our proposed scheme 
is superior to (about 2-dB) the detector in IfTOl . The main reason is that the differential detection 
approach employed in iflOll was based on the estimation of the previous symbol. Consequently, 
when one symbol was decoded incorrectly, it will affect the decoding of the consecutive symbols 
thus leading to serious error propagation. Comparatively, the information about the estimation 
of the previous symbol is not required in our proposed differential detection and is, thus, able 
to prevent the error propagation. 

In Fig. [3l we include the Genie-aided results by assuming that each source node can perfectly 
remove its own information from the received signal. It can be noted from the results that the 
proposed differential detection scheme introduces negligible performance loss compared to the 
genie-aided scheme. We also compare the BLER performance of the differential scheme over a 
GSM channel (a practical channel) and a quasi-static Rayleigh fading channel. From the figure, 
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it can be observed that there is almost no performance loss in a GSM channel compared to 
the quasi-static Rayleigh fading channel which clearly justifies the robustness of the proposed 
differential scheme in slow fading channels. It also indicates that the effect of non-constant 
channel on proposed scheme can be ignored which validate our assumption of quasi-static fading 
channel model. 

In Fig. BJ we show the optimum power allocation scheme of the DDSTC-ANC scheme. It can 
be seen that more power should be allocated to Pi when the channels from relay nodes to T 2 
are better than the channels from relay nodes to T x . It is interesting to find that when a 2 ^ a 2 , 
the sources still share half the total power for the optimal power allocation. 

In Fig. [51 we examine the BLER performance of the proposed scheme with power allocation 
for the system with four relay nodes. The SORC is used at relays and signal is modulated from 
a BPSK constellation. We also take into account the relay's location as: case 1 (the symmetric 
case), where relays are placed halfway between the source nodes, i.e., Ti,T 2 , and a 2 = 1 and 
a 2 = 1; and case 2 (the asymmetric case), where relays are close to the source node T 2 , and 
a 2 f = 1 and a 2 , = 10. It can be observed from Fig. [5] that the BLER performance of the proposed 
scheme with power allocation can provide considerable performance gain in comparison with 
the equal power allocation (EPA) scheme, i.e., Pi = P 2 = Pr. = 

VII. Conclusion 

In this paper, we have proposed a DDSTC-ANC scheme for TWRNs with multiple relays. A 
simple differential signal detector was developed to recover the desired signal at each source. The 
performance of the proposed DDSTC-ANC scheme was analyzed and the OPA was presented 
to improve the system performance. Analytical results have been verified through Monte-Carlo 
simulations. 
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Fig. 1. Block diagram of the two-hop TWRN. 
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Fig. 2. Simulated BLER performance using Alamouti codes (2 relays). 
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Fig. 3. Simulated BLER performance using Alamouti codes (2 relays) over a GSM channel and quasi-static Rayleigh 
fading channel. 
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Fig. 4. Optimum power allocation between source and relay nodes. 
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Fig. 5. Simulated BLER performance by the proposed DDSTC-ANC using SORC with transmit power allocation (4 
relays). 



